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ABSTRACT 


(iii) 


This work preserits analytical and ntiaierical models for the 
analyaivS of flow and bed level vaiiation in gradual open channel 
expansions with movable beds. The analytical solutions are 
obtained for the equllibt'ium bed profile and corr espond i ng flow 
parameters based on the one dimensional flow assumption. Cartesian 
coordinate as well avS a cylindrical coordinate system (radial 
flow) are used for this purpose. The equilibrium bed profile 
correspond.s to steady flow and no furt.her aggrada t loti . Tlie 
numerical model solves the governing partial differential 
equations for- water flow and sediment trans|>ort by 

finite-di f f ererice methods. A quasi-steady uncoupled approach is 
adopted, whetein tht^ flow et|aations and the e( 4 uations for sediment 
transport are solved separately. Flow conditions are assumed 
steady during the solution of equations for sediment flow. 
One-dimensional as well as two-dimensional models are developed. 
Tlie twO“d imens i onal models solve for flow parameters by the false 
transient method . Two numerical schemes , Lax-diffusive and 
MacCormack vUcdieme are used . A vSimple algebraic coordinate 
transformat i on technique is used to convert the expanding physical 
domain into a rectanguiat- compu t a 1 1 onal domaiii to facilitate easy 
application of boundary conditions. The numerical models predict 
both temporal and spatial variation in bed levels. All the models 
presented here are valid only for uniform &lze sediment and for 
flows without separation. A comparison l^etween the numerical 
results and the analytical results for the equilibrium bed profile 
shows that the numerical models correctly solve the governing 
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equations. A s ens i t i vi t y analysis is alvSO performed to study the 
numerical effects. A parametric study is conducted to investigate 
the effect of various parameters on aggradation in open channel 
expansi ons . 
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CHAPTER I 


Introduction 


An open channel transition is a channel reach in which 
cross-section or slope changes with distance such that there is a 
change in the flow parameters (Henderson 1966). Channel expansions 
and contractions are common examples of such transitions. They are 
usually required in human made channels for several practical 
purpos es . They are used in structures between canals and 
f lumes , canals and tunnels and canals and inverted syphons to 
reduce energy losses(Chow 1959). Many hydraulic structures such as 
barrages and bridges may require constriction of the channel in 
order to reduce the length of the structure and thereby the cost. 
These structures require contractions on the upstream side and 
expansions on the downstream side. This is especially true in case 
of sediment-laden flows because channel bed may aggrade or 
degrade depending on the flow conditions. Understanding of flow 
phenomenon in these transitions and ability to quantify the 
variations in flow parameters is essential to the proper design 
and operation of hydraulic structures such as syphons, aquaducts, 
barrages, bridges etc. Expansions and contractions may also occur 
naturally. For example, rivers expand and widen before they enter 
reservoirs, oceans or other main rivers. They also widen when they 
leave steep mountainous areas and enter gentle plains. Flow 
analysis in sucli natural expansions helps in our understanding of 
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reservoir sedimentation and delta formation. 

In an expansion, the flow area increases in the down- 
stream direction and as a consequence the velocity 
deer eases . Therefore , the sediment carrying capacity gets reduced as 
the water enters an expansion and this results in sediment 
depos i t ion . An expansion built with a uniform slope initially, will 
aggrade till it attains an equilibrium bed level. Such a state la 
reached when the sediment carrying capacity is same at every 
croaa-section in the expans i on . The above simple picture gets 
complicated for the case of large angle and sudden expansions . The 
flow separation and the attendant recirculation zones change the 
flow characteristics and sediment movement signi f icant ly . Inf act , 
Nashta et al,(1987) have observed erosion at the entrance region 
of a sudden channel expansion with movable bed. The process 
gets even more complicated if the sediment is not uniformly 

sized. The larger particles tend to move towards the wallCNashta et 

al.l987) and the finer particles get deposited downstream, During 
degradation, finer particles are transported in preference to 

larger particles. There may be a stage when the bed is 

predominantly made up of large particles which may provide a 
protective layer for the underlying finer particles. This process 
is called armouring of bed. 

The analysis of flow and bed variation in expansions is 

complicated because the flow is two-dimensional . Difficulties are 

/ 

also encountered if the flow is supercritical and shocks are 
present. The governing partial differential equations do not have 
analytical solutions for majority of cases and therefore, require 
numerical solution. Application of boundary conditions in these 
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numerical schemes is not straight forward and may pose 
difficulties. The time scale of flow variation is significantly 
different from the time scale of bed level variation and this 
should be taken into account while developing an efficient 
numeriacal algorithm. Even if the numerical scheme is efficient 
and reliable, uncertainty in the computational results arises 
because the relationship between water and sediment still defies 
rfgorous treatment from the point of view of basic mechanics. 

The present work is an attempt to study the bed level variations 
in open channel expanslotis with movable beds. It is limited to 
subcritical flow in gradual expansions where there is no flow 
separation ahd to beds comprising of uniform size sediments. 

1. 1 REVIEW OF LITERATURE 

In the past, several attempts have been made to study the 
flow characteristics in channel expans ions . Chat urvedi ( 1 9 63 ) and 
Garde et al.(1979) have studied the subcritical flow in a gradual 
expansion , Lokrou and Shen (1983) analyzed the flow in sudden 
expansions by a similarity approach while HehtaC1979) studied the 
flow characteristics in two-dimensional expansions. Hinds (1928), 
nitra(1940) and Chaturvedi (1963) proposed different design 
procedures for expansions with given flow conditions and certain 
assumptions based on experiments .Vittal and Chiranjeevi (1983) 
proposed a rational method of design for open channel 
transitions -ChowC 1959 ) ,Scogerboe et al.C1971) and Garde et 
al.(1979) attempted to design devices to control separation in 
expanding channels . 
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Host of the studies on open channel expansions are with 
respect to fixed-bed channels . Studi es on movable bed ew pan s i ons are 
only few. Ashida and Hiyai(1964) experimentally • analysed the 
formation of deltas in an abrupt expansion .Bates (1953) studied 
the bar formations in streams as they come out of a natural 
constriction ChangC1982) and Hert ens ( 1 986 ) ( Raudkivi 

1990) described experimental data on aggradation in the form of 
delta formation in reservoirs. It is found that deltas in 
reservoirs grow both in length and width. The initial development 
is a strongly elongated delta in flow dir ect ions , over which the 
flow depth shallows. This encourages sidewards flow. The other 
effect is rise in the water level in the channel upstream of the 
reservoir. The reservoir in the above experiments is modelled as a 
sudden expansion. Figure 1.1 shows the formation of sand delta in 
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laboratory experiments . Recent ly , NashtaC 1984 ) and Nashta et 

al-(1987) investigated subcritical flow in abrupt expansions with 
movable beds. They found that the high velocity jet scours the 
material of the bed. The curvature of the separating stream line 
makes the coar-ser sediment to move and deposit towards the walls. 
The finer material is transported downstream and is deposited in 
the form of a low transeverse bar. Figure 1.2 shows the flow 
pattern and the bed topography in a sudden expansion with movable 
bed. The available experimental work in movable bed expansions is 
mostly related to sudden expansions and studies on gradual 
expansions are not many. 

Laursen (1958,1963) derived analytical solutions for scour depth 
in long contractions. He solved the flow equations as given by the 
Manning's formula along with a sediment discharge relationship for 
the condition of final equilibrium bed profile. His solution is 
valid only for predicting the maximum scour depth. 

It does not describe the general two-dimensional bed level 
variation. Similar approach may be used for channel expansions. 
Straub(1953) also used similar approach but with a di f ferent 
sediment transport relation. 

Previous studies of open channel expansions have been 
mostly either experimental or dealt with the analytical solutions 
to simple cases. Also, no attempts have been made to study the 
temporal variations of bed level in an expansion duo to 
non-uniform flow conditions .Although experimental and analytical 
studies contribute significantly to our understanding of open 
cha;nnel flows with movable beds, mathematical modelling (using 
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numerical ahemes) of movable bed hydraulics has become an active 
area of research in computational hydraulics .Hathematical models 
for movable bed channels essentially solve the governinfi equations 
for water and sediment flow. Water flow is usually represented by 
the ah^low water equations which describe the conservation of 
mass and momentum. The sediment flow is represented by the sediment 
continuity equation and a sediment transport equation which 
relates the sediment discharge to the flow parameters . This 
sediment transport equation is generally an empirical equation 
based on laboratory and field investigations . The above governing 
equations for water and sediment flow constitute a set of 
non-linear hyperbolic partial differential equations for which 
analytical solutions are possible for only simple 
cases . 

Ther ef or e , they are solved using numerical schemes . Either 

f init e— di f f er ence or finite-element methods can be used for this 
purpos e . Hos t of the standard movable bed models can be classified 
into either • uncoupl ed models, wherein the water-flow equations and 
sediment continuity equation are uncoupled during a given time step 
or coupled models, wherein all the governing equations are solved 
s imul tan eous ly . They can also be classified as either unsteady 
models or quasi-steady flow models. The unsteady models solve the 
complete Saint Venant equations describing the water flow. On the 
other hand, quasi— steady flow models assume water flow to be 
steady during the computation of bed level variation .Models also 
differ with regards to the physical conception and the numerical 
scheme employed for the solution. Some models use uniform size 
sediment while others use non-uniform size sediment and attempt to 
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simulate armouring process also. For the modes of transport of 
0 ediment , some use total load concept while others make a 

distinction between different mode of transport “• A detailed 

classif icat ion of models is presented in figure 1.3 . 


MinthermatiesiX mod^X 



Figure 1.3 Model classification 

Although the achievements In the fixed bed 
modelling (bed level is atfssumed invariant) are considerable , same 
can not be expected of movable bed models. The fixed”"bed open 
channel flow is well described by partial differential equations 
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developed by Saint Venant in 1871. The numerical solution of these 
equations and the results obtained therefrom have consolidated the 
base of computations in fixed-bed flows. On the contrary, the 
physics of the movable bed flows is poorly understood .The 
understanding is not satisfactory even for conditions like 
one-dimensional flow with uniform size sediments. In practical 
situat ions , the sediment particle size is non-uniform and the flow 
is at least two-dimensional if not three-dimensional . Besides , the 
armouring and sorting effects may be present . However , considering 
the time scale, the cost and difficulty in measuring the flow 
parameters and bed levels in field investigations .mathematical 
modelling is sometimes preferred to experimental 

invest igat ion . Also , mathematical models are more convenient to use 

and several alternatives can be studied easily. But, the model has 

m 

to be r el labl e , robust , accurat e and economical (Holly 1986). 

Since the early work of Vreugdenhil and de Vries 
in 1967 (Cunge et al . 1980) , several one-dimensional movable bed 
models have been developed , Hoi lyCl 986 ) and Dawdy and Vanoni Cl 986 ) 
presented excellent reviews of numerical simulation of alluvial 
hydraulics , HEC-6 .developed by U.S. Army Corps of 

Engineers, FLUVVIAL-SCChang and Hill 1976) , lALLUVIALCKarim and 
Kennedy 1982) and CHAR-2 (SOGREAH 1978) are some of the 

commercially available movable bed models for one-dimensional 
flows.HEC-6 is a good example of a quasi-steady flow model whereas 
CHAR-1 (Cunge et al . 1980) is an unsteady model. 

One-dimensional movable bed models are 
many . However , only few two-dimensional movable bed models are 
available. The models developed by de Vriend and StruiksmaC 1 983 ) 
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a.nd Shimizu and ItakuraC 1989 ) are applicable only for curved 
channels .Ilatsutomi et al.(1984) developed an unsteady model for 
routing dam break floods in natural channels with movable 
beds . However , i t has not been applied to study the long term bed 
level changes in gradually expanded channels . TABS- 2 (Waterways 
Experiment Station 1987) is a finite element model for water and 
sediment routing .It solves unsteady shallow water equations and 
sediment continuity, and is primarily designed for analysing 
sedimentation in large bodies of water like estuaries .reservoirs 
etc. It is not efficient for analysing long term bed evolution in 
channel transitions which are usually not very long. 

l.a SCOPE OF THE PRESENT STUDY 

The literature review presented in the previous section 
shows that majority of experimental studies on channel transitions 
with movable beds are for sudden expansions . There are only few 
studies on gradual expansions . Although a number of mathematical 
models are available for movable bed hydraulics , most of them are 
for one-dimensional flows. The available two-dimensional models are 
either not applicable to the case of movable bed expansions or 
inefficient. 

The scope of the present work is to 

1) develop- a two-dimensional mathematical model for analysing 
long term bed evolution in gradual expansions 

2 ) develop an appropriate one-dimensional movable-bed model which 
can be applied to gradual expansions 

3) develop analytical solutions to the equilibrium bed profile in 
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an expansion by assumina the flow to be radial so that the 
mathematical models can be validated and 

4) to study the aggradation characteristics in a gradual expansion 
through a parametric study. 

The numerical models presented in this study use a 
quasi-steady uncoupled approach for simulating long term bed»^ 
evolution in gradual expansions . The false transient method is 
employed to obtain the steady flow parameters . Two numerical 
schemes, Lax-diffusive scheme and HacCormack scheme are used for 
this purpose. The application of the models is limited to 
subcritical flows in gradual expansions where there is no flow 
separation. Also, the bed is assumed to be made of uniform size 
particles . 

Chapter II presents the governing equations and the 
assumptions made in deriving them. Analytical solutions are 
derived for few simple cases in chapter III . The numerical 
schemes for the solution of the governing equations are described 
in chapter IV. Chapter V discusses the computational results and 
the conclusions are presented in chapter VI . 



CHAPTER II 


Governing Equations 


To simulate the unsteady flow in open channels with 
movable bed, the governing partial differential equations for the 
flow of water and sediment are numerically solved. Although 
extensive research has been carried out to understand the exact 
relationship between water flow and sediment movement C Shen 
1970, Garde and Rangaraju 1985), the present knowledge in this area 
can only be considered semi-empirical. Also, complete solution of 
three-dimensional equations of motion for water and sediment is 
very complicated. Since the aim of the present study is to 
simulate only the bed level variation in a channel expansion, the 
following two major assumptions are made in the analysis. 

(1) Ue are not interested in computing the velocity 
of sediment particles and (2)the sediment discharge at any point 
is uniquely related to the depth averaged flow parameters at that 
point . 

The above assumptions make the sediment continuity equation 
sufficient to describe the sediment flow .The equation for the 
conservation of momentum for the sediment is implicitly 
represented by the relationship between the sediment discharge and 
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depth averaged flow parameters. 


In this c h a p t e r r g o v e rn i n q 
sediment are presented. Co-ordinat 
transforming the physical domain of an 
rectangular computat i onal domain is 
equations for one-dimensional assumption 


equations for water and 
e transf ormation for 

expanded channel into a 
discussed. SiiTiplified 
are also discussed. 


2.1 WATER FLOW 


The two-dimensional unsteady gradually varied flow 
equations in open channels are (Lai 1977, Jimenez 1987) : 

continuity equation 

Sh . <9(Vx h) . «9CVy h) _ „ (2.1) 

momentum equatioii in x-direction 

, ■»( V;; vy h ? , .3 

St « K ' 

( 2 , 2 ) 

momenttim equation in y-direction 

£1^) ^ '>q« , yy ... „ , hj . ^ h , (,h^ , 2 } , t so - sf ) 

at ax *»y * ^ y y ^ 

(2.3) 


in which, h is the flow depth, Vx and Vy are the depth averaged 
velocities in x-direction and y-direction respectively ,g is the 
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acceleration due to gravity, and SO and SO are bed slopes in x 

K y 

and y- directions respect ively. x and y are co-ordinate axes, and t 
is time. 

The friction slopes are calculated using the following equations , 


SF = 

K 


2 


n 


Vk 


ivx + 

. 1- 333 
h 



C2.4) 


SF = 
y 


2 

n 


Vy 


i V« »• Vy 

, 1- 333 
h 


C2.5) 


in which, n ia Hanning’s roughness coefficient. 


F 

Although n is in general a comlicated function of flow 
depth, bottom roughness , slope , discharge and bed forms, a constant 
value is assumed throughout this study. However elaborate 
equations can be easily incorporated because explicit numerical 
methods are used for the solution of governing equations. 


Equation 2.1 - 2.3 are obtained by depth averaging the 

thr 6 e“dimens i onal equations and making the following assumptions 

( 1 )acceleration in the vertical direction is 
negligible. In otherwords pressure distribution along the * 

v<a.Tr-Kc<s>Jl i® hydrostatic. 

(2)the velocity distribution is uniform over the flow 
depth. 

C3)bottom shear stress is dominant and all other shear 



15 


stresses are negligible. 

C 4 ) the friction losses computed using steady state 
formulae are valid even for unsteady flow conditions. 
C5)the channel bottom slope is small. 

The above assumptions are valid for most of the gradually 
varied flow si tuations . However , the governing equations do not 
account for the effective stresses which arise due to 
Ci)lamlnar viscous stresses, 

(ii ) turbulent stresses and 

( i i i ) stresses due to depth averaging 

Extra turbulent .stress-like terms appear while depth 
averaging the momentum equations because of the non-uniformity of 
the velocity in the vertical direction .Based on experiments in 
the laboratory channels, Odgaard and Bergs(1988) have shown that 
the error introduced by uniform velocity assumption is negligble. 
Flokstra(1977 ) also showed that away from the walls the effective 
stresses are dominated by the bottom stress . However , it should . be 
noted that these effective stresses should be considered while 
simulating circulating flows (Flokstra 1977). Therefore, models 
presented in this study are not valid for large angle expansion 
where flow separation occurs. Few models are available for 
simulating effective stresses using turbulence closure 
models . However , these methods, as applied in open channels , (Rastogi 
and Rodi 1978) are at best successful! for fixed bed 
channels. Also, their application may be required only if one is 
interested in three-dimensional flow structure and actual sediment 



16 


movement . 


2.2 SEDIMENT FLOW 


As dlacuased earlier, only sediment continuity equation 
and a sediment discharge-water flow relationship are reqiiJteci for 
completely representing the sediment flow. In cartesian co-ordinate 
system the sediment continuity equation is given by 


^ . _1_ 
St 1-X 




^y 


= 0 


( 2.63 


Where, Z is the bed elevation, gs^ is sediment discharge 

per unit length in x-direction , gs^, is sediment discharge per 

unit length in y-direction and A ia the porosity of the 

se<t i m e n t par t i o 1 - e e . g a and gs depend on the flow parameters at 

M y 

that point. Several relations are available for estimating these. 
These relations are well documented elsewhere (Graf 1972, Garde and 
Rangara Ju , 1 985 3 and are not described here. As suggested by 
Vanoni(19753 and Jansen et al. (19793 many of these equations can 
be represented in the following functional form. 


X = f ( Y 3 


(2.7) 


Where ,X is a transport parameter and Y is a flow 
parameter * The transport parameter depends on the sediment 
discharge and the grain properties .while the flow parameters 
depend on the flow properties and the grain diameter. The numerical 
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models in the present study are structured in such a way that any 
sediment discharge relationship can be easily considered. However, 
the following equation is adopted for its simplicity. 






az/ 2 


C2.8) 


Uhere, ga is the resultant sediment discharge per 
unit length, and a.^ , a^are empirical constants which depend on 
sediment grain properties, bed forms and flow properties. Although 

p 

the above relationship is a highly simlified one , it has been 
used in mathematical models by many researchers to gain insight 
into the problem of aggradation and degradation (Gill 1983, Zhang 
and Kahawital987 , Park and Jain 1987 ) .Analytical solution for bed 
level variatiohs in a channel expansion could be obtained by the 
use of above equation. This made it possible to verify the 
numerical models by comparing the numerical results with the 
analytical solution. Before closing this discussion, it should be 
mentioned that the above equation is valid only for uniform 
sediment size. Therefore, the models do not give proper results in 
case of degradation where armouring effect is present. During 
degradation, finer particles are transpo’tsed in preference to 
larger particles. There may be a stage when the bed has only 
sufficiently large particles which can provide a protective layer 
for the underlying finer particles. This protective layer is 
refered to as the armoured layer. It should also be mentioned at 
this stage that no distinction is made between the bed load and 
the suspended load and equation 2.8 is assumed to represent the 
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total sediment discharge. The aim of the study is to develops a 
numerical model for analysing aggradation in channel expansions 
and validate the numerical model by comparing the numerical 
results with the anlytical results. Therefore, the use of a 
simplified relationship is justified. However, a more suitable 
sediment discharge relationship should be included and the model 
should be calibrated before it is used in field situations. As 
mentioned earlier, any sediment discharge relationship can be 
easily incorporated into the models presented in this study. 

a. 3 CO-ORDINATE TRANSFORMATION 

Equations 2.1, 2.3 and 2.6 are for a cartesian co-ordinate 
system. Therefore, it becomes necessary to replace the boundaries 
of an expansion by an approximate grid (figure 2.1) to solve the 
governing equations by finite difference techniques. 


Ac tual 
boundary 



Approxirr'iaic 
boundaf y 

r 


Flow 



FIGURE 2-1 BOUNDARY APPROXIMATION 
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Figure 2.2 Coordinate transformation 


Thia approximation may adversely affect the solution (Roache 1972), 
This may be taken care of by transforming the physical domain to a 
rectangular computational domain (figure 2.2) using the following 
algebraic transformations ( Anderson et al. 1984) 

(2.9) 

( 2 . 10 ) 
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Where, b(:<) is the distance between the symmet ry 1 ine and the 
upper boundary (figure 2 . 2 -a). 0 nly half of the transition is to be 
corriputed because of symmetry . Now the boundaries of the 
computational domain coincide with 17 = 0 and T) = 1 (figure 
2.2-b).The governing equations need to be transformed so that tliey 
can be used with reference to the rectangular computational 
plane .Applying the chain rule of differentiation 



or 




>7 b a 

b ) ‘ »r) 


( 2 . 11 ) 


where, b C? ) = - 4- ^ 


Sirailarly, 

& _ 1 ^ 

■" bCf ) <»r) .12 ) 


Substitution of equation 2.11 and 2.12 in equations 2.1, 2.2 ,2.3 

and 2.6 and subsequent simplifications lead to 



continuity equation for water 


[ t>(? ) h] & [ b(^ ) Vx h ] <» fVy h-T) b C^D Vx hl-O 

dt &K ^ ^ * 

(2.13) 

momentum equation for water in x-direction 

hiK ) (h Vx^+ fth ?2)J + 

^ [vx Vy h - Y),b(?) (h Vx^+ eh ^2)J 


^ [b(?) Vx h 1 S 

3^ * ^ 


t 


= g h • [- {Z b(^)} +1^ {17 b' (?) Z) - SF-^] 

(2.14) 

momentum equation for water in "J-direction 

£_[bC?) Vy h ] . £_r b(?) Vx Vy hi 

St d? I J 

+ [Ch Vy2+ g.h^/2) - Vx Vy h T) b (? ) ] = b(? )g SfJ 

(2.15) 


continuity equation for sediment 


( 1-X )bC?) 


sz 

3t 




- « b « ) ■ e».] . 0 

(2.16) 





Note that the bed slopes SO and SO are replaced by the 

^ y 

following 

C2.17D 

on - 

SOy - ~ C2.18) 

It should also be noted that the physical domain is not a function 
of time i.e. the water surface width is not changing with time. 
Therefore, the application of the model is limited to channel 
cross-sections with vertical side walls. 

2.4. ONE-DIMENSIONAL GOVERNING EQUATIONS 

The one-dimena ional governing equations can be obtained 
from the two-dimensional equations by averaging them for any 
cross-section and assuming that the velocity, the bed level and 
flow depth do not vary across the width. These are given as 
below: 


continuity equation for water 


3ft 


+ 




0 


C2.19) 


momentum equation for water 


^ 0^/A ) 

^t dx 


+eA^=gA(SO-SF) 


C2. 20) 



sediment continuity equation 
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In the above equations , Q = flow discharge, G = 
sediment discharge and B = channel width. These are well suited 
for rectangular channels with high values of inlet to outlet width 
ratio. However, the assumption is not valid for lower values of 
inlet to outlet width ratio. In this case a better approximation 
would be to assume that the flow is radial. Refering to figure 2.3 



V 


FiGUf^E 2-3 f'F f !f II T Of- (RADIAL 

FLOV/ 
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the one-dixnena ional radial flow equations are (Abbot 197 9 , Townson 
and Salihl 1989) 


Continuity equation for water 


dh 


Sr 


+ 


. Vr h 
Sc 


0 


(2.24) 


momentum equation in r-direction 

SVr „ SVr Sh SZ ^ OCX 

Wr 'S'r ^ ° (2.25) 

continuity equation for sediment 

= 0 

sediment discharge equation 

G = ( Vr )^2 (2.27) 

the resistance relation will be 

n2 2 

SF = (2.28) 

2^2^3.333 
r p n 


where, Vr is the velocity in radial direction, r is the radial 
coordinate and (9 is the angle of expansion. Though the 
one-dimensional equations are strictly valid for uniform flow 
condition across the width or across the arc at a radial 
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distance, as the case may be, they still provide satisfactory 
results for flow and bed variation in channel transitions. Also, 
the equilibrium or steady bed level variation in the channel 
transition can be obtained by integrating the steady 
one-dimens ional equations. These analytical solutions of the 
steady one-dimensional equations are presented in the next 
chapter - 



CHAPTER III 


Analytical Solution 


The one-dimensional flow and sediment equations are amenable to 
analytical solution for simple boundary and steady state 
conditions. A channel expansion with an initially uniform bed 
slope as shown in figure 3.1 is considered . Even if the downstream 
and upstream flow conditions i.e. the downstream flow depth and 
the flow discharge as well as the sediment inflow into the 
channel are not changed , the channel bed will evolve with time 
due to the presence of non-uniform flow conditions. The 

velocity and flow areas are varying in the downstream 
direction. As a consequence , the sediment inflow into any 
control volume is not same as the sediment outflow. This imbalance 
results in either aggradation or degraddation of the channel bed 
as the case may be . These bed level changes continue till the 
channel bed attains an equilibrium state such that the sediment 
carrying capacities of all the sections along the channel are 
equal. The steady one-dimensional flow and sediment equations can 
easily be Integrated to obtain the above final equilibrium bed 
profile in an expansion. The derivations based on radial flow as 
well as rectangular flow assumptions are presented in the 
following sections. 




28 


3.1 RADIAL FLOW IN AN EXPANSION 


For steady conditions equations 2.24-2.26 become 


dV» 

dr 


+ V>' 


dh 

dr 


= 0 


C3-1) 


Yr 


dVr 

dr 


+ 



dZ 
* dr 


+ ft SF 


0 


C3.2) 


dG 

dr 



(3.3) 


Equation 3.1 can be written as 



h 

r 


or 


d [ Vr h] 
Vr h 



(3.4) 


Integration of equation 3.4 leads to 

log { Vr h > + loft r = loft C, 

IP* ^ ,, !!*«*«••••••••••••*••********* ^ 

Where. is the constant of integration. 
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sine e , 


0 = Vr h C r ) ; 


Vr h r 


(3.6) 


similarly, equation 3.3 gives 


G r = C, 


(3.7) 


where, is an integration constant, 

Substitution of equation 2.27 in equation 3.7 leads to 


ajVr r = 


(3.8) 


h r 


'2 

r = C, 


(3.9) 


h = C 3 r 


l/a^ - 1 


(3.10) 


a 1/a 

where, = [ g-j - ] “ 77 “ 


(3.11) 


at the downstream end,r = r^ and h = h^ where, h^is flow depth 
at the downstream end. 

Substituting the above boundary condition in equation 3.10, 


we get 


^l^a„-l 


(3.12) 
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equation 3.2 may be written as 

dZ _ _gp _ ^ _ Vr dVr 
dr dr fi dr 

from equation 2.28 


-SF = 

where 


2^2. 3- 333 
r ft h. 


3- 333 

1- 333- 

V ^2 


= - 




equation 3.10 gives 

dr ^7 r ^ 

where Cy ■= - C^Cl/a 2 “^^ 

from equation 3.15 

_ Vr ^ -2^a -1 

g dr 8 ^ ^ 

2 


.(3.17) 

. (3,18) 

.(3.19) 

(3.20) 

(3.21) 

(3.22) 



31 


where Cg 




C 


aubstitution of equationa 3.18,3.20 and 3.22 in 
leada to 


dZ 

dr 




1- 333- 


3. 333 


+ C_ + C- r" ^ V ^ 

f o 


equat ioa 


(3 


Int^gratiag equation 3.24 with respect to r » we get 


. 2. 333-3. 333^* y-1. -2^^ „ 

Z C5,r + C^gr + C^^r + 


C3 


•' Where, 

^9 “ 2.333-3.333r 

C 7 

^10 ~ r -1 

c 

11 2 r 

r = l/a^ 

and Cj ^2 integration constant. 

At r “ (r^ is the radius of the inflow section), Z “ ZI 

the bed elevation at the inflow section). 

Applying the above boundary condition to equation 3.25, we 


C 

C 

C 


CZI 

get 


2 


^ 2. 333-3- 333y 

ZI + CpT + 

r 2 . 333-3. 333r . 

“ L 



r?: 


-1 


'10" i 



C 


(3 


3.23) 

3.17 

. 24) 

.25) 

) . 26 ) 

1.27) 

1.28) 

is 

.30) 


Equations 3.13,3.15 and 3.30 describe the steady state variation 
of flow depth , velocity and bed elevation in an expanding channel 
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with movable bed. These are based on one-dimensional radial 
assumption. The equations for rectangular flow assumption 
presented next. 

3.2 RECTANGULAR FLOW IN AN EXPANSION 


For steady state conditions, equations 2.19,2.20 
2.21 are simplified to 


dx 


4^2-^ A 


dx 

dx 


dh 

dx 


A (SO - SF) 


[b a^(0/A)^2] = 0 


(3 

C3 

(3 


Integration of equation 3.31 gives 

Q = S 3 

where is a constant of integration. 


A “2 

or BC-a2)A~^2"^ ^ + a"*2. B = 0 


where, B == 


dB 

dx 


or 


-Ba^(l/A)Cl/A^2)[ ||B + hB]+A^2.B = 0 


(3 


From equation 3.33 

a, (Q)*2 

1 ^ ^ dx , a„ ■* : 


C3. 

.(3 


or dh 
dx 


(r-i)B 


Ik. 

B 


dCQ^/A) . , dh _ _ 

^ ^ ^ = « 




2.23 in 

equation 3. 

r_ 15 _ _ 

^2 2 1 

On I 

1. ci X 

i:v J 


C3 


flow 

are 


and 

.31 ) 
.32 ) 

.33 ) 

.34) 

.35) 

.36) 
.37 ) 

.38) 


f 3 .39 'I 
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d Z 

where, ~ dx (3.40) 

or, ::iQ^ [^B + hB']+^ = -4^ qIbI— (3.41) 

gA3 ^ dx ^2^3. 333 • 

or, dZ _r Bq3 1 ^ Q^h B _ 

dx “L 3 ^ J dx ,3 ^2^3.333 (3.42) 

gA gA B h 

Subatitution of equation 3.38 in equation 3.42 leads to 
dZ _ Cy-l)B'h j. BO^ Q^hB* O^n^ 

dx - B L .3 ^ J ^ ~3 ^273.333 (3.43) 

gA gA B h 

Integration of equation 3.38 gives 

h = B^"^ (3.44) 

where, is an integration constant. 

Application of the boundary condition h =h^at B = BO ( BO = width of 

the downstream end ) results in 




Integration of equation 3.43 is not straightforward if not 
impossible. However, it can easily be integrated numerically by 
using either Hodified Euler’s method or Fourth order Rangekutta 
method. Uhile doing so the boundary condition Z = ZI at B = BI is 
incorporated. Modified Euler's method and Fourth order Rangekutta 
method for steady open channel flow conditions are explained in 
detail by Subramanya ( 1986 ) and Chaudhry ( 1990 ) and therefore, 
they are not repeated here. 



3.3 CLOSURE 


3 4 

In thia chapter analytical solutions for steady state 
equilibrium bed profiles in channel expansions have been 
presented. However, analytical solutions for the time variation 
of bed levels are not possible and numerical methods need to be 
applied . The numerical schemes for one dimensional as well as two 
dimensional cases will be discussed in the next chapter. 



CHAPTER IV 

Numerical Solution 


The governing equations for unsteady flow of water and 
sediment transport as discussed in chapter II constitute a set of 
non-linear hyperbolic partial differential equations ( Lyn 1987 ) . 

Analytical solutions for these equations are available only 
idealised cases . Therefore, they are solved numerically .In this 
chapter , numerical schemes are presented for the solution of 
one-dimensional as well as two-dimensional bed level changes with 
respect to time . 

I 

Bed level changes occur at a much slower speed 
than the flow changes and unsteady terms in water flow equations 
may not be important while simulating the bed transients ( de Vries 
1965 ) . This is especially true for the cases studied here where 

change in bed level occurs due to non-uniformity of flow parameters 
rather than the unsteady variation of velocity and depth 
Therefore , a quasi-steady uncoupled model can be used for 
predicting the long term bed level changes which occur in a channel 
expansion . The quasi-steady state means the flow parameters for 
water are constant during a small period of bed elevation change 
due to sediment flow , This model is more efficient than the 
unsteady coupled model from the point of view of computational 
costa . In a quasi-steady , uncoupled model , the equations for 
water flow and sediment transport are solved separately . The 
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strategy for solution ( Fi^ 4.1 ) is to first compute the steady 
flow parameters for computing the sediment discharge values at 
different sections . These sediment discharge values are then used 
in sediment continuity* equations for computing the bed level 
changes . After solving the sediment continuity equation , the flow 
computation is repeated with the new bed configuration and the 
procedure is continued untill an equilibrium bed profile is 
reached. . 

4.1 COMPUTATIONS FOR ONE-DIMENSIONAL MODEL 


For steady state conditions ^ equation 2.20 may be 
modified and written as 


7 * 9 7 

^ _ Olb ^ Qlji 

dh dx .3 „2 , 3.333 

_ &k B h 

dx Q^B ■_ - 


( 4 . 1 ) 


Since Z,Q and B values are known at any x , the right 
hand side is a function of only h and x . The above equation can be 
integrated by Ilodified Euler’s method to obtain the h values at 
any cross section . Subsequently , the velocity in any cross 
section can also be determined . The numerical integration starts 
at the downstream end since the h value is known at this section 
and then proceeds in the upstream direction . The above computed h 
and V values can be used to determine sediment discharge in each 
cross section using the equation 2.22 . The sediment continuity 




Floiv chart for qiiasi-slcady uncoupled model 


F 
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equation ( equation 2.21 ) is represented in the finite difference 
form to compute the bed level changes as given below . 


* 



C4 .2) 


(4.3) 


where , 


F = B 




(4.4) 


In equations 4.2 , 4.3 and 4.4 , the subscript i refers 
to the node number in the finite difference grid and the 
superscript * refers to the time level t+At . Substitution of 
equations 4.2 , 4.3 and 4.4 in equation 2.2 and subsequent 
simplification results in the following : 


^2 

At a^O 

■ 


Ax B^Cl-X) 

h®2 

L 1 

J 


(4.5) 


Equation 4.5 can be used to determine the bed level at 
time t+At at all nodes for i=2,3..,N+l where N is the number of 
reaches into which the channel is divided . The bed level at node 
1 la given by the inflow boundary conditions . In the present study 
a straight reach is included in the upstream of the expansion 
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and since there is no change in the sediment inflow , the bed level 
at node 1 does not vary with time . The new bed levels at time 
t+At are substituted in equations 4.1 and the procedure is repeated 
to obtain the bed levels at time t+2At . The computations are 
continued till there is no appreciable change in the bed level and 
the channel attains an equilibrium state . 

4,. 2 COMPUTATIONS FOR TWO-DIMENSIONAL MODELS 

The two dimensional equations are solved using the same 
principle as for the one dimensional equations . However , the 
steady flow parameters in two dimensional model are obtained using 
the false transient method as described by Bhallamudi and 
ChaudhryC 1991 ) . In this method , the unsteady water flow 
equations are solved to obtain the steady state values using time 
as an iterative parameter . Any initial conditions through the 
channel are assumed and the boundary conditions are set equal to the 
actual steady state conditions . Using the two dimensional unsteady 
flow equations , the system conditions are computed for a 
sufficient length until the variation of flow variables are ■ 
negligible and the conditions have converged to the steady state 
values corresponding to the given boundary conditions . Note that 
the time in this method is used as an iterating parameter and 
therefore , the computational time step here is not same as the 
computational time step At used in the sediment continuity 

equation. In the present analysis , two different finite 

difference methods are used for water flow computations . These 
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J 


(4.9) 


In equations 4.6 to 4.9 F represents any 
variable of which the partial differentiation is approximated . In 
equation 4,9 , F represents the source terms as averaged by the 
values at neighbouring nodes . The subscripts l,j refer to the grid 
points in the ? m\d r) directiotis respectively . The superscripts k 
and k*^l refer to the values of the variable at time level k and k+1 
respectively . k is the known time level while k-*-! is the unknown 
time level . Substituting equations 4.6 4.9 in equtaions 2.13 

2.15 , algebraize equations are obtained and are solved 

explicitly for the flow paramaters at the unknown time level . 


4. 2.1.1 Initial Conditions 


To start the unsteady state calculations 

V , V - h and Z at time t = 0 are to be specified 

nodal points . Since a false transient method is used 

ayteady state conditions , approximate values of V , V 

X y 

given as Initial conditions . 


values of 
at all the 
to obtain 
and h are 
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4- 1 . (S Bouridary conditions 

There are three types of boundaries to be considered 
in this analysis . These are i) open boundaries , ii) symmetrical 
boundary and iii) solid side wall boundary . 

Open boundaries 

An open boundary could be either an inflow or an outflow 

boundary . The specification of boundary conditions depends on the 

nature of the flow . Since only sub-critical flows are considered 

in the present study, two boundary conditions need to be specified 

at the upstream end and one boundary condition at the downstream 

end . Flow depth at the down stream end and discharge or at the 

upstream end constitute these boundary conditions . A straight 

portion is included upstream of the channel and the flow is assumed 

to be one-dimensional at the upstream end , Therefore , the 

transverse velocity , is zero at this point and this forms the 

second boundary condition , Even after the specification of 

boundary conditions from the point of view of well posedness of the 

problem , extra equations need to be used to completely define the 

values at the upstream and downstream ends . This is due to the 

finite difference requirements . Strictly speaking , characteristic 

method ( Chaudhry 1987 ) should be used to determine V and V at 

X y 

downstream end - However, a simple extrapolation from the values 
already calculated at the interior nodes using Lax-diffusive Scheme 
gave satisfactory results .Such an extrapolation procedure has been 
used successfully by Bhallamudi and Chaudhry C 1991 ) . 
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Symmetric boundary 
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Reflection procedure is adopted for 
boundary ( Roache 1972 ) .Refering to 
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Figure 4.3 Reflection technique for symmetric boundary 
imaginary reflection points is considered .All non-conservative 
flow variables other than the normal velocity are specified as even 
functions with respect to the symmetric line while the normal 
velocity is specified a-s an odd function . This makes the average 
normal velocity at the boundary equals to zero . 

Solid side wall boundary 


For the models presented in this study, the slip 
condition is the proper boundary condition at a side wall since 
only bottom shear stess is considered . Therefore , the resultant 
velocity at a solid wall is tangent to it . Anderson et al.(1984), 
have discussed several boundary techniques in gas dynamics 
applications. The reflection principle seems to be most suitable 



for the present applications. It should be remembered that the 
reflection technique is only approximate for solid side walls 
Refering figure 4.4, the flow depth and the magnitude of the 



Figure 4.-4. Reflection technique for solid side wall boundary 
imaginary velocity at reflection point are same as that of 
corresponding interior grid point . In order to make the normal 
velocity at wall equal to zero , the velocity components at 
reflection point are computed using the following formulae . 


Vxp= VcoaCZe - «) (4.10) 

Vyp= VsinCZ^ - «) (4.11) 

where , 

V= Resultant velocity at the interior point , 
x-component of the velocity at the reflection point , angle 


between the wall and x-direction and « = angle between the resultant 
velocity V and x-directlon . 



4. a. 1.3 stability 
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Lax-di f f ualve Scheme haa to eatiafy the Courant - Friedrichs - Levy 
(CFL) condition for stability . The above conditions in case of 
transformed equations with transformed co-ordinates may be 

expressed as : 


'2d 


AtC V + ■/ fih ) 
b ( X ) Ay) 


. a^2 + 


{ b(x) Ay) } < 1.0 C4 .12 ) 


equation 4.12 has to be satisfied at every grid point for the 
scheme to be stable. The above condition is heuristic and is based 
on a linearised form of the governing equations for one-dimensional 

i 

flow. It also does not consider the effect of friction and bed slope 
terms . Therefore , a certain amount of numerical experimentation is 
required before choosing the upper limit of . 

A.Z.Z Setdlmenl flow > 


As already mentioned ,the algorithm uses a 
quasi-steady flow assumption and therefore, the equation for 
sediment continuity is solved separately after obtaining the steady 
flow parameters . The following finite difference approximations are 
used to solve the sediment continuity equation . 
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where, F is any variable whose partial derivative needs to be 
determined. Substitution of equations 4.13 - 4.15 in equation 2.16 
gives an explicijt equation for 'Z’ at the unknown time level. 

It should be noted again that the At used inequation 4.13 
is not same as At used in the false transient method for water 
flow. The computational time step used for solving the sediment 
continuity equation is much larger than the At used in the 
Lax-diffusive Scheme. There is no theoretical work regarding the 
stability of the quasi-steady uncoupled models for bed-level 
changes and the computational time step in equation 4.13 is chosen 
after an extensive numerical experimentation.lt should be as large 
as possible from the point of view of computational efficiency, and 
at the same time should lead to stable and convergent results. 

Before closing this section, it should be mentioned 
that the upstream boundary condition for the sediment continuity 
equation is same as that for one-dimensional model . Bed level at 
the downstream end is obtained by extrapolation from the interior 


points . 



4 .. 2. 3 MacCoriaack sch<£rme 
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In the present study, the water flow equations 
are solved by HacCormack Scheme also. This is done to ascertain the 
accuracy of numerical results. HacCormack Scheme is a second order 
accurate explicit scheme. It comprises of a predictor part and a 
corrector part. In the predictor part, forward finite differences 
are used while in the corrector part, backward finite differences 
are used for approximating the spatial differential terms . Ref ering 
figure 4.2 the predictor and corrector parts are defined as 
below . 

fBfedictor part f 
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Corrector part 
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Uhere, F is any variable and i and j refer to 
nodes in ? and 7 ) directions respectively. Superscript^ k is for 
known time level, ^ is for predicted values and ** is for corrected 
values. The value of 'F' at the unknown time level k-*-! is given by 
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Several alternatives for MacCormack scheme are 
also possible. One can use backward finite difference in corrector 
part for calculating the spatial variabl es . Another alternative is 
to change the direction of differencing from one time step to the 
next. The initial conditions, boundary conditions and stability 
criteria are same as discussed for Lax-diffusive Scheme. 

In this chapter the numerical schemes were described for 
the solution of bed level changes and water flow parameters with 
respect to time. Computer programmes are prepared using the above 
'principles. These programmes are used for analysing sediment 
deposition in channel expansions- The results of the study are 
presented in the next chapter. 



CHAPTER V 


Results and DtscussioNS 


In the present study, a two-dimensional, quasi-steady 
uncoupled model has been used to analyse the variation of flow 
parameters and bed levels in a gradual expansion with movable 
bed. As discussed earlier, most of the experimental works in 
gradual expansions are for fixed beds while most of the studies in 
movable beds are for sudden expansiojia . Also , moat of the work 
pertains to the study of equilibrium bed profiles. Experimental 
data (Son! et al . 1980, Brush and Uolman and — B e p - l-and 1960, Begin 
et al . 1981) for temporal variation of bed levels due to unsteady 
boundary conditions are available only for one-dimensional 
case . Experimental data are not available for bed level variation 
in gradual expansions. Therefore, the accuracy of the numerical 
models developed in this study is checked by comparing the 
numerical results with the analytical results for one-dimensional 
case derived in chapter II I. Once again it should be borne in mind 
that proper implementation of boundary conditions and selection of 
the time steps used in flow equations and sediment continuity 
equation are important aspects for successful mathematical 
modelling . 

In the present study, five models have been used. For 
brevity in discussion, these models are refered to as 
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(I)IIODEL 1 : Analytical solution of one-dimensional equations in 

cylindrical coordinates for equilibrium bed profiles. 

(ii) MODEL 2 : Analytical solution of one-dimensional equations in 

cartesian co-ordinates for equilibrium bed profiles. 

(iii) MODEL 3 : Numerical solution of one-dimensional equations in 

cartesian coordinates for quasi-steady flow and unsteady bed 
conditions . 

(iv) MODEL 4 : Numerical solution of two-dimensional equations 

for quasi-steady flow and unsteady bed conditions using 
Lax-diffusive Scheme for solving water flow equations. 

(v) MODEL 5 : Numerical solution of two-dimensional equations for 

quasi-steady flow and unsteady bed conditions using 
MacCormack Scheme for solving water flow equations. 

Models 1,2 and 3 are one-dimensional while models 4 and 5 
are two-dimensional. Models 1 and 2 give analytical solution and 
models 3,4 and 5 give numerical solutions . Therefore, models 2 and 3 
form the link between analytical and numerical models. 

In the following sections , model verification, sensitivity 

J 

of the numerical models to numerical parameters and a parametric 
study for final equilibrium bed profile are presented. 

5.1 MODEL VERIFICATION 

Model verification has been done in two phases. First, 
all the one-dimensional models are compared among themselves 
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(b) C ross-sec i ion 

FIGURE 5-1 CHANNEL EXPANSION 


(models 1, 2 and 3) and then, the two-dimensional numerical results 
are compared with the one-dimenalonal results (models 2,4 and 5). 





5.1.1 Input Values 
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The analytical as well as numerical results have been 
obtained for the following input values . Figure 5.1 shows the plan 
and cross-section of the channel considered. In this figure, 

Lj= ength of inlet portion = 5.6 m, 

Lg= length of expansion = 4.0 m, 

Lq= length of outlet portion = 6.4 m, 

BI = channel width at the inlet, 

BO = channel width at the outlet = 0.9 m and 

BCx) = channel width at any x. 

Although the numerical models can be applied to any type 

of expansion, only a straight walled expansion is considered 

here. The inlet channel width is varied from 0.5 to 0.8 m in the 

numerical exper iments . The initial uniform bed slope in the 

dfpwnstream direction, S0^= 0.001 while the transverse slope, SO = 

0.0. The downstream control depth, h.= 0.107 m and flow discharge , 

d 

Q= 0.019 m /s -The sediment transport parameters, a^= 0.000145 .and 
a^^B.O . The bed level at the upstream end, 21 = 0.0 and the 
porosity of channel bed, ^ = 0.4 . The above geometric, flow and 
sediment properties roughly correspond to the experimental set up 
of Nashta(1984) for sudden expansion with movable bed. 

The numerical models used a finite difference grid 
of A^ = 0.4 m and Ayj = 0 - 1 - A Courant number value of 0.85 is used 
to achieve numerical stability in the false transient method of 



solving water flow equations. The computational time step, At for 
the solution of sediment continuity equation is 1000 seconds. 
However, At has been varied upto 5000 seconds to study its effect 
on numerical stability and accuracy. 


5. 1 . 2 Orie-dliwerislonal Models 

Figures 5.2 to 5.4 show the variation of flow 
depth, bed level and deposition obtained by using one-dimensional 
models for BI = 0.8 m 1 26 ) and BI = 0.6 m C/^= 4°18 ).The 
flow depths for models 1 and 2 Cfifture 5.21 are matching while 
model 3 underpredicts depths for small ft values . However , the 
predictions are satisfactory for higher ft values or lower BI 
values (figure 5.2 b). Results from models 1 and 2 are matching 
very well because of small value of ft, in which case the radial 
and rectangular flow assumptions are equivalent. The error in 
predictioiis by model 3 may be due to the numerical approximations. 
However, the error is not significant and ammounts to less than 
one percent . 


The variations in the bed levels along the channel 
for the final equilibrium state are shown in figure 5.3. It can be 
observed that the numerical results (model 3) compare 
satisfactorily with the analytical results (model 1 and 2). As 
before, the comparislon is better for larger values of ft (=» 4 ) 

than for smaller values of ft (=k1'“). The numerical results vary in 

I 

a stair case fashion. Similar type of variation can be observed 



FLOW DEPTH (m) FLOW DEPTH (m) 




(■ a j ■ 1 



DISTANCE (m) 

MODEL 1 + MODEL 2 O MODEL 3 


Chi ft m 4 *^ 18 * 



DISTANCE (m) 

MODEL 1 <- MODEL 2 O MODEL 3 


Figure 5.2 Comparison of one -dimensional models for flow depth 


variation 
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even in flood routins by Lax-diffusive scheme. This . is a 
characteristic of the finite-difference scheme adopted. It can 
also be observed that the rectangular flow model (model 2) 
overpredicts the aggradation as compared to the radial flow model 
(model 1) for large values of ft (figure 5.4 b) . However, this 
difference is not very significant . It should be noted that the 
results for models 1 and 2 are presented along the channel center 
line. The actual bed level contours as obtained using radial flow 
model will differ slightly from those obtained using the 
rectangular flow model. This is illustrated in figure 5.5. 

The above discussion shows that the concept of 
quasi-steady uncoupled modelling as incorporated in model 3 can be 
successfully applied to predict the long term bed level variation 
in gradual expansions. However, model 3 assumes one-dimensional 
flow conditions and should not be used in situations where flow 
is highly two-dimensional. In the next section, results obtained 
using two-dimensional models are presented. 


^ - - *. 

1 2 , 5 



rectangular flow — ■ radial flow 


Figure 5.5 Bed level contours 



5-1.3 Two-dlwiensilorial model ss 


The two-dimensional models (model 4 and model 5) also use 
the concepts of quasi-steady flow conditions and uncoupling of 
sediment and flow equations As mentioned earlier, they use the 
false transient method to obtain the steady flow parameters 
during any computat ional time . Model 4 uses Lax-diffusive scheme 
and model 5 uses MacCormack scheme for this purpose. Two different 
numerical schemes are used in this study so that there is more 
confidence in the numerical results . Figure 5.6 - 5.8 show the 
results obtained- using these models. The analytical results 
obtained using model 2 are also presented in these figures for the 
purpose of comparislon. As before, all the computations are made 

o' o' 

for two different expansion angles, ft - 1 26 and ft - 4 18 . The 

comparislons are made for variation along the centerline as well 
as along the wall for the final equilibrium state. 

Figure 5.6 shows the variation of flow depth with 
distance . As can be observed from these figures, the numerical 
results from models 4 and 5 compare very well with each other. The 
numerical results also compare satisfactorily with the analytical 
results obtained using model 2. Although flow depths obtained by 
one-dimens ional model are generally higher than those predicted by 
two-dimensional models (figure 5.6), this difference is not 
significant. Figure 5.7 shows the variation of bed levels with 
distance for the final equilibrium state. As before, the numerical 
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results obtained uain^ model 4 match well with those obtained 
using model 5 giving certain amount of confidence in our numerical 
results. The stair case type variation can be observed in the case 
of model 4 which uses Lax-diffusive Scheme. This is eliminated in 
case of model 5 which uses MacCormack scheme It should be noted 
that MacCormack Scheme is second order accurate while 
Lax-diffusive scheme is only first order accurate. The 
two-dimensional effects in aggradation can be observed from figure 
5.7 . Considering the deposition along the wall, the bed levels 
obtained using models 4 and 5 are higher than the bed levels 
obtained using the one-dimensional model in the initial reaches of 
the expansion. As a consequence, the two-dlmens ional models give 
lower bed levels than those obtained by one-dimensional model 
towards the outlet of the expansion. However, this difference is 
not significant along the center line, especially for the smaller 
ft values (figure 5.7 a-ii). It can also be observed that, as 
expected, the two-dimensional effects become stronger when the 
expansion angle ift)i& large . Figure 5.8 shows the variation of 
deposition with distance. Similar trends are observed as before. 
Although the numerical results obtained using the two-dimensional 
models differ from those obtained using one-dimensional model, 
they are of the same order of magnitude. This illustrates the 
satisfactory performance of two-dimensional numerical models. 

B. 2 SENSITIVITY ANALYSIS 

The developement of any mathematical model is not 
complete untill' a sensitivity analysis is performed for dependence 
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on numerical parameters. In this section, the sensiti:vity of 
numerical results to i^rid size and computational time step are 
discussed. The results are presented for the same flow and sediment 
data as before- Although computations are done for both 

. o. ' o' 

and for = 4 18 , results for only =1 X8 are presented. 
Similar results are obtained for ft = 4^16 .The numerical results 
obtained using model 4 as well as model 5 are considered. 

5.^.1 Effect of Grid Size 

Figure 5.9 shows the effect of computational 
distance step ftn downstream direction, DS on the deposition for 
model 4. As can be seen ,the difference in the numerical results 
for the three values of DS = 0. 2,0.4 and 0.5 is not significant, 
especially for deposition along the center line. However, this 
insensitivity is not as good for deposition along the wall. This 
may be due to the reflection boundary technique used for 
simulating the side wall boundary conditions. This technique is 
only approximate and unfortunately, no other alternative method is 
available .Similar observations can be made for numerical results 
obtained using model 5 (figure 5.10). Figures 5.11 and 5.12 show 
the sensitivity of numerical results to the value of computational 
distance step in transeverse direction, DN for models 4 and 5 
respectively. Three different values of DN equal to 0.833 , 0.1 
and 0.125 are used in the computations . Again, same a-a-m -e trends 


as discussed above are observed. 
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S. 2. 2 Effect of Computational time step 

No theoretical guidelines are available for 
selection of the computational time step, DTINE in quasi-steady , 
uncoupled models. Therefore, it is necessary to do some numerical 
experimentation before a proper value is chosen. It should be 
chosen such that it is not too large so as to introduce 
instabilities . At the same time , it should not be too small so 
that the computational cost is not prohibitive. Figures 5.13 and 
5.14 show the effect of computational time step on numerical 
results for models 4 and 5 respectively. Three different values of 
DTIME equal to 1000 s, 2000 s and 5000 s are used. As can be 
observed from these figures, all the three DTIHE values gave the 
same results. Although not presented here, an extensive numerical 
experimentation indicated that the maximum value of DTIME that can 
be used without stability problems strongly depends on the 
expansion angle as well as sediment parameters. For example , 

O ' ^ 

stable results could not be obtained using model 5 for = 4 jlS 

when DTIME was greater than 2000 s . The results presented here 
indicate that once we get stable results with a particular DTIME 
value, any DTIME value smaller than this would not increase the 
accuracy significantly. 


5.3 PARAMETRIC STUDY 

In this section, the effect of different parameters 
on the final equilibrium state is studied. Only one parameter 
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is vai'ied at a time and others are kept constant. Table 5.1 
the several cases studied 


shows 


SI. Mo. 

BI 

TABLE 5.1 

^1 

^2 

^d 

1 

CO 

o 

0.000145 

5 . 0 

0.107 

2 

' 0.7 

0.000145 

5.0 

0.107 

3 

0.6 

0.000145 

5.0 

0.107 

4 

0.6 

0.000290 

5 . 0 

0.107 

5 

0.6 

0.000435 

5.0 

0.107 

6 

0.6 

0.000145 

5.5 

0.107 

7 

0.6 

0.000145 

4 . 5 

0.107 

8 

0.6 

0.000145 

4 . 0 

0.107 

9 

0.6 

0.000145 

5.0 

0.117 

10 

0.6 

0 . 000145 

5 . 0 

0.127 


5.3.1 Effect of sediment transport parameter o-j 

Referinfi, to equation 3.30 , the final equilibrium 
levels in the expansion do not depend upon a.^. The same can 
observed from fieure 5.15 . Here , three different values of 
equal to 0.000145, 0.000290 and 0.000435 are used. 


bed 

be 

The 


equilibrium bed profile is same for all the three cases. 
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(o) ALONG WALL 



DISTANCE (m) 


(b) ALONG CENTER UNE 



Figure 5.15 Effect of on bed elevation 
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B« 3. 2: Ef'fect of iSediirfeGrrit, Transport. Param&for a^ 

Figure 5.16 shows the effect of on final 

equilibrium bed levels along the wall (figure 5.16 a) as well as 
along the center line (figure 5.16 b). Four different values of 
equal to 5.5, 5.0, 4.5 and 4.0 are considered. It is clearly 
indicated that larger the value of a^ smaller the value of 
deposition or aggradation. If the flow parameters are such that 
the velocity is less than one, then a higher value of a« indicates 
less movement of sediment. Consequently, the aggradation will be 
less if all the other conditions are same. In the present case, 
the flow velocity is of the order of 0.2 m/s ~ 0.1 xn/s . This 

explains the trends observed in figure 5.16 . However, the effect 
is highly non-linear as expected. 

I 

5. 3« 3 Effect of downstream flow depths 

Figure 5.17 presents the variations of final equilibrium 

bed levels along the wall and along the centerline for three 

different cases of h. equal to 0.107 m ,0.117 m and 0.127 m . For 

d 

the computations made here, the downstream depth does not seem to 
have significant effect on tlie equilibrium bed profile. However, a 
higher downstream depth seems to increase the deposition (figure 
5.17 b). A higher depth means a lower velocity and conseqently 
increased deposition. 










BED LiVEL <m 
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r 7 9 

6.3.4. Effect of expansion ratio 

1 

Three different Inlet widths, BI equal to 0.8 in, 0.7 m 
and 0.6 in are considered. Figure 5.18 shows the effect of BI or 
expansion ratio on the aggradation. As expected, higher expansion 
angle or a lower BI value results in moce deposition and 
consequently higher equilibrium bed levels. As can be observed 
from figure 5.18, the deposition increases in the downstream 
direction . 

6.4. TEMPORAL VARIATION OF BED LEVELS 

The models presented in this study are verified by 
comparing numerical results with analytical results for final 
equilibrium bed levels. However, the main advantage of the 
numerical models is their ability to predict the temporal 
variation of bed levels. It is expected that the models give 
proper results even during the unsteady state. Unfortunately, it 
is not possible to verify this due to lack of analytical results. 
For the sake of completeness the temporal variation of channel bed 
profile is presented in figure 5.19 . Bed levels at time , t equal 
to 0.0 seconds (initial condition), 5 days, 25 days , 50 days and 
100 days are shown. The final equilibrium bed profile is obtained 


at t = 100 days . 



o 
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Figur* 5.18 Effect of expansion ratio on bed elevation 
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Figure 5.19 Temporal change in bed elevation 
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CHAPTER VI 

Summary atid Recommendations for Furtfier Investigation 


6.1 SUMMARY 

This investigation presented analytical and numerical 
models for the analysis of flow and bed level variation in gradual 
open channel expansions with movable beds. Analytical solutions 
were developed for the equilibrium bed profile for steady 
conditions ahd no further aggradation. On the other hand, 
numerical models can predict both temporal and spatial variation 
of bed levels. A qu<tjBi~8teady , uncoupled approach was adopted in 
developing the one-dimensional as well as two-dimensional 
numerical models. The one-dimensional numerical results matched 
very well with the analytical results which are also for 
one-dimensional case. Two different numerical schemes were used in 
developing the two-dimensional models. Results from these two 
models matched with, each other very well. The results were also 
comparable to those obtained using one-dimensional analytical 
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model. This indicated the satisfactory performance of the models 
as far as the solution of igovernirn^ equations was conducted. 
However, the #^iimple sediment transport equation used in the 
models forms a weak link. The sediment transport parameters should 
be properly evaluated before the models are applied in field 
situations. The models presented in this study were structured in 
such a way that any empirical sediment transport equation can be 
easily incorporated. A parametric study usin« the numerical models 
indicated the following. 


1.) The sediment transport parameter have no erffect on the 
equilibrium bed profile. 

2) An increase in the parameter decreases the equilibrium bed 
levels if the flow velocity is less than one. 

3) The downstream flow depth slightly increases the deposition. 

4) As expected, larger the expansion angle, greater is the 
deposition. 

6.2 RECOMMENDATIONS FOR FURTHER WORK 

The following is a list of subjects that require further 

research. 

1) A completely unsteady coupled model needs to be developed for 
analysiiig the aggradation and degradation of bed levels in 
transitions during flood flows. 
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2) Several different boundary conditions should be included in the 
model for analysis of flow and bed in any ijeneral transition and 
for any upstream conditions like sediment over-loading. 

3) A proper method for including the effect of armouring due to 
non-uniform size sediment needs to be explored. 

4) There is no experimental data regarding the aggradation in 
channel expansions and there is an immediate need for such work. 



REFERENCES 


1. Abbot, n.B.,( 1 97 9 ), Conipu t a t i onal Hydraulics; Elements of the 
Theory of Free Surface Flow, Pitman Publishing Limited, 
London . 


2. Anderson, D, , Tannehill,J. and Fletcher, R.,(i9843» 
Computational Flu.id Hechanics and Heat Transfer, 
Hemisphere Fublisliing Co. 


3. Ashida.K. and Hiyai , H . , ( 1 9 6 4 ) , "St udy on the Sedimentation in 
an Abrupt Expansion Annual Report of DPRI, Research 
Institute for Disaster Prevention, U.O kyoto, Kyoto-fu, 
Japan . 

4« Bates , C . C . , ( 1 9 5 3 3 , "Rational Thpory of Delta Formation," 
Bulletin of the American Association of Petroleum 
Geologist , Vol . 89 ,no . 9 , pp .2119-2162. 


5. Begin, Z.B., Heyer, D.F. and Schumm, S.A., (19813, "Development 

of Longitudinal Profiles of Alluvial Ctiannels in 
Response to Base -Level Lowering." Earth Surface 
Processes and Land Forms, vol. 6, no. 1, pp . 49-60. 

6. Bhallamudi, H . S . and Chaudhry, ft . H . , ( 1 9 9 1 3 » '’'Computat ion of 

Flows in Open Cliannel Transitions", Jour. Hydr . 
Res.,IAHR,to be published 


,7. Brush, L.n. and Uolman, M.G., (19603, "Knickpoint Behaviour in 
Noncohesive Platerial: A Laboratory Study." Bulletin of 
the Geological Society of America, vol. 71, pp. 59-74, 


0. Chang, H.H.,( 1982 }, "Hathernat leal Hodel for Erodible Channels," 
Jour. Hydr. D i v . , ASCE , vol . 108, no. 5, pp . 678-688. 


9. Chang, H.H., and Hill, J.C., (19763, "Computer modelling of 
erodible flood chaniiels and deltas," Jour. Hydr. 
Div.,ASCE, vol. 102, no. 10, pp. 1461-1477. 


10. Chaturvedi ,ri.C., ( 1 9 6 3 3 , F low Charact erst ics of Axi-symmetric 

Expans ioti ASCE, Vol. 8 9,Hyd-3, pp . 61-92. 


11. Chaudhry ,11 . H .,( I 987 ) , Applied Hydraulic Trans i ents , 2nd ed. , 
Van Nostrand Re inhold Company, New York. 



12 . 


Pr ent ice 


Hall , 


Chaudhry , tl , H . , ( 1 9 9 0 ) , Op en Charm el FI ow , 

Eng 1 ewood , M . J . 


13- Chow, V.T., ( 1959), Open Channel Hydraulics, ricGraw-Hill Book 

Co . , New York . 


14. Dawdy, D.R. and Vanoni, V. A., (1986), ^'Kodellinfi Alluvial 
Channels,’' Water Resources Research, vol . 22, no. 9, pp . 
71&-813. 


15. de Vriend, H.J. and Struiksma, N.,(1983), "Flow and 

Deformation in River Bends , "River ' 83, ASCE, New 

Orleans, October 24-26, pp . 991- 1011 


16. de Vr i ea , H. ( 1 9 6 5 ) , "Considerat ions about ^^on-Steady Bed-Load 
Transport in Open Channels," 11 Confiress, lAHR, 
Leningrad , Paper 38 


17. Flokstra, C., (1977), "The Closure Problem for depth, averagrd 

Two-Dimensional Flow,” 17 Congress lAHR, Baden-Baden, 
paper A106, vol. 2, pp . 247- 256. 


18. Garde, R . J . , Rangara ju , K.G. anrJ Hishra, R.C., ( 1 979 ), 

"Subcritical Flow in an Open Cliannel Expansion," Jour. 
Irr. and bower, Central Board of Irr*. and Power, 

India, Vol. 36, no. 1 

\ 


19. Garde, R.J., and Rangara ju , K . G .,( 1 985 ) , Mechanics of Sediment 
Transportation and Alluvial Stream Problems, 2nd. ed., 
Wiley Ea.stern Ltd., New Delhi, India. 


20. Gill, n. A., (1983), "Diffusion Model for Aggrading Channels,” 
Jour. Hydr. Res. vol. 21, no. 5, pp . 355- 367. 


21. Graf, U.H.,(1971), Hydraulics of Sediment Transport, McGraw 
Hill Book Co. , New York- 


22. Henderson, F. n.,(1966), Open Channel Flow, MacMillan Co., New 
York. 


23. Hinds ,J.,(1928), "The Hydraulic Design of Flume and 

Transitions," Tran.s., ASCK, Vol, 92, pp . 1423-1485. 

24. Holly Jr., F.M.,(1986), "Numerical Simulation in Alluvaial 

Hydraulics,” 5th Congress of the Asian and Pacific 
Regional Division of the lAHR, August 18- 21, Seoul, 

Kore 



25.Janaen, P.Ph., Bendefiom, V.L., Berg, V.D.J., de Vries, H. and 
Zaneu, a. (1979), Principles of River Engineering, The 
Non-Tidal River, Pitrnan Publishing Ltd., London, 


26. Jimenez, O.F., ( 1987 ), "Computation of Supercr i t i cal flow in 

Open Channels," thesis submitted for the degree of 
riaster of Science, Department of Civil and 
Environmental Engineering, IJashington State University, 

27. Karim, H.F. and Kennedy , J . F ,,( 1 982 ) , "lALLUVIAL" :A Computer 

Based Flow-and-Sediment-Rout ing Hodel for Alluvial 
stream , an Application to Hissouri River, lAKR 
Report No. 292, University of Iowa City. 


20. Lai, C,,C1977)," Computer Simulation of Two-dimensional 
Unsteady Flows in Estuaries and Embayments by the method 
of Charac t er i s t i CvS , U.R.I. 77“-85,U.S. Geological Survey 


29. Laursen, E. N . , ( 1 9 5 8 ) , "Th e Total sediment Load of Streams, 
Proc. ASCE,vol. 108, HYll, 1282- 1295 . 


30. Laursen, E, N . , ( 1 963 ) , "Analys ivS of Relief Bridge Scour", Proc. 
ASCE, vol.84,HYl, 1-36 


31. Lyn, D. A., ( 1987 ), "Unsteady Sediment Transport tlodel ing" , Jour , 
of Hydr . Engr . , ASCE ,Vol. 113, no. 1, January, pp . 1-15. 


32. Lokrou, U.P. and Shen, K.U., ( 1 983 ) , " Anal ys is of the 

Char act ers t i cs of Flow in sudden Expansion by Similarity 
Approach," Jour, of Hydr. Research, lAHR, Vol . 21, no. 

2, pp, 119-131. 


33. Hatautomi, H,, Asada, H. and Sato, T.,(1984), 
"Two-dimensionalHodeling of Dam Break Floods in Natural 
Channels with a Novable Bed", Channels and channel 
control s true t ur es , ( Ed . K.V.H. Smith), Springer 
Verleg, pp . 4.89- 4.101. 


34. !1ehta,P.R. , (1979) , "Flow Chare t ers tics in Two-Dimensional 

Expansion," Proc. Jour, of Hydr. Div., ASCE, vol. 105, 
no, Hy-5, pp. 501-516. 


35. Ultra, A. C., ( 1 9 4 0 ) , Techni cal Hemorandum No. 9, UPRI(India). 


36.Nashta, C.F,,(i984) Subcritical Expansions in Open Channels 
with Hovabl.e beds, Ph.D. Thesis, Univ, of Roorkee. 



37. Mashta, Garde, R.J. and Swamee, P . K . , ( 1 987 ) , '"Subcr i t i cal 

Flow in Open Channel Expansions with Movable Bed”, 
Jour. Hydr. Res., vol, 25, no. 1, pp. 89-102. 


38. Odfiaard, J.A. and Ber£s, M. A., (1988), "Flow Processes in a 
Curved Alluvial Channel", Uater Resources Research, 
vol. 24, no.l pp . 45-56. 


39. Park, I. and Jain, S. C.,C1986), "River Bed Profiles with 
Imposed Sediment Load", Jour. Hydr. Engr . , ASCE, 
vol. 112, no. 4, pp. 267-279. 


40. Raatofii, A. K. and Rodi, U . , ( 1 9 78 ) , "Pr edi c t i one of Heat and 
Mass Transfer in Open Channels" , Jour . Hydr. En^r., ASCE, 
vol. 104, no. 3, pp. 397-419. 


41.Raudkivi, A . J . , ( 1 9 9 0 ) , Loos e Boundary Hydraulics , 3rd 

ed.,Perfiamon Press . 


42, Roache, P. J.,(1972), Computational Fluid Dynamics, Hermosa 
Publ ishers . 


43, Sco^erbo e , G . V . , Austin,L,H. and Bennel.R.S., (1971), "Energy 
loss Analysis for Open Channel Expansion," Jour, of 
Hydr. Div,, Froc. ASCE, Vol. 97, Hy-10, pp . 1 7 1 9 - 1 7 3 6 . 


44. Shen, H. U . ( Ed . ) , ( 1 9 7 1 ) , River Mechanics, vol‘,V , Water 
Resources Publication, Fort Collins, Colorado. 


45. Shimizu, Y, and Itakui'a, T.,( 1989),” Calculationa of Bed 
Variation in Alluvial Channels", Jour. Hydr, Eofir., 

ASCE, vol, 115, no, 3, pp, 367-384. 


46. SOGREAH Consulting En^i neers , ( 1 97 8 ) , Mod el in^s of River 

Longitudinal Profile Evolution Usinfi the CHAR-2 System. 


47. Soni, J,P., Garde, R.J. and Raju, K.G. R . , ( 1 980 ) , " Afigradat i on 
in Streams due to Over loadiri^s” » Jour. Hydr. Div., ASCE, 
VO 1 . 1 0 6 , no. 1 , pp . 117-132. 


48. Subramanya, K.,(19S6), Flow in Open Channels , Tata 
McGraw-Hill Publishing Co. Ltd., New Delhi, INDIA 


49. Straub, L. riisnouri Rivhi Report, U.S. A51imy Corps 

of Ensineeir., House document 2 38, Appendix XV, pp. 1156. 



50. Townson, J.Il. and Salihi, A. H . , ( I 989) , "riodels of Dain-Break 
Flow in R"T Space'" , Jour . Hydr . Engr.^ASCE, vol.ll5j no. 
5, pp. 561-575 . 


51, Vanoni , V, A., ( 1975), Sediinentat ion Encineering, ASCE- Hanuals 
and Report, 9 on Engineering Practice, no. 54. 


52. VittaljN. and Chi ran j e evi , V . V , (1983), "Open Channel 

Transitions: Rational Method of Design," Jour, of Hydr. 
Div., Proc. ASCE, Vol . 109, NO. 1, 


53. Zhang, H, ahd Kahawita, R .,( 1987 ), "Non-linear Model for 
Aggradation in Alluvial Channels", Jour. Hydr. Engr . , 
ASCE, vol. 103 , no. 3 , pp . 



